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1 Introduction 



Quantum field theories with a topological term ( "#-term" ) in the action have proved 
to be particularly challenging to investigate. Such theories are related to a few 
important open problems in theoretical physics, including the so-called "strong CP 
problem" in strong interactions, and to interesting phenomena in condensed matter 
physics, such as the quantum Hall effect (for a recent review on theories with #-term, 
see Ref. pQ). 

On the one hand, topological properties are intrinsically nonperturbative, thus 
requiring a nonperturbative approach to the study of these systems. On the other 
hand, the most effective of these approaches, namely the numerical study by means 
of simulations in lattice field theory, cannot be directly applied to these systems, 
due to the presence of a so-called sign problem. Infact, the complex nature of their 
Euclidean action prevents the computation of the relevant functional integrals by 
means of the usual importance-sampling techniques. Numerical investigations have 
then required the use of techniques which allow to avoid the sign problem, usually 
based on analytic continuation or on the resummation of the contributions of the 
various topological sectors to the partition function [21 |3j HJ EJ EJ U\ E] • The basic idea 
of these techniques is to modify or split the functional integral, in such a way that 
the resulting expression(s) have a positive-definite integration measure, and therefore 
can be treated with the usual numerical techniques. The difficulty of dealing with 
an oscillatory integrand is however not completely overcome, but simply shifted to 
the problem of reconstructing the original functional integral, which is usually a 
very delicate issue from the numerical point of view. It is worth noting that, beside 
having their own theoretical interest, theories with a #-term share the sign problem 
with finite-density QCD, and so the development of techniques and algorithms to 
solve or by-pass the sign problem can have positive consequences on the study of 
the QCD phase diagram by means of numerical simulations. 

Among the various existing models, the two-dimensional 0(3) nonlinear sigma 
model with #-term (0(3)eNLcrM) deserves particular interest. It has been shown 
long ago by Haldane |9] [10] that chains of quantum spins with antiferromagnetic 
interactions, in the semiclassical limit of large but finite spin S, are related to this 
model at coupling g 2 = A/[S(S + 1)], and at 9 = or tt if the spin is respectively 
integer or half-integer. Haldane conjectured that quantum spin chains for half- 
integer spins show a gapless spectrum, and correspondingly that a second-order 
phase transition takes place in the 0(3)gNLaM at 9 = it, with vanishing of the 
mass gap and recovery of parity. Arguments supporting this conjecture have been 
provided in Ref. [TT]. Moreover, in Ref. JT5] it has been argued that the critical 
theory for generic half-integer spin antiferromagnets is the SU(2) Wess-Zumino- 
Novikov-Witten (WZNW) model [131 US EES] at topological coupling k — 1, which 
in turn should determine the behaviour of the mass gap near 9 = n. 

Numerical investigations of Haldane's conjecture have been performed, following 
basically three different strategies. A first strategy [21 [TS] is based on the determina- 
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tion of the probability distribution of the topological charge by means of simulations 
at 9 = 0, which allows in principle to reconstruct the expectation values of the var- 
ious observables at 9 ^ 0. In order to achieve the very high accuracy required by 
this approach, the authors of Refs. [21 [16] have employed a constrained ( "topolog- 
ical" [17]) action on a triangular lattice, which allows simulations by means of an 
efficient Wolff cluster algorithm [18] . The parameters of the action were chosen in 
order to be in the weak-coupling regime. Using finite size scaling theory, the authors 
of Ref. [21 EE] found a second order phase transition at 9 = ir, in agreement with Hal- 
dane's conjecture, and a finite size scaling in good agreement with the assumption 
of a WZNW-type of critical behaviour. 

A second strategy [19] is based on the determination of the mass gap at imaginary 
values of 9, that can be obtained directly by means of numerical simulations, and 
the subsequent analytic continuation to real values of 9, in order to check if the mass 
gap vanishes at some point. The authors of Ref. [19] found indeed that the mass 
gap vanishes at 9 = 9 C for some real 9 C , and moreover that 9 C = ir within the errors, 
again in agreement with Haldane's conjecture. 

Finally, the third strategy [201 [2B [221 [23] makes use again of numerical simula- 
tions at imaginary values of 9, in order to determine the topological charge density, 
and of a controlled way of performing the analytic continuation to real 9 that greatly 
reduces the uncertainties connected to this process. Applying this strategy to the 
CP 1 model, that is expected to be equivalent to the 0(3) model, the authors of 
Ref. [22] found a richer phase structure, with a first-order phase transition at 9 = tt 
for < 0.5, and a line of second-order phase transitions with recovery of parity for 
0-5 J$ J$ 1-5, with continuously varying critical exponent. At ~ 1.5 the critical 
exponent becomes 2, and parity is recovered analytically. 

In this paper we want to investigate further on this issue, by applying the strategy 
of Refs. [201 [2H [221 123] directly to the 0(3) e NLaM. Our aim is to understand the 
origin of the discrepancy between the results of Refs. [21 ES] and those of Ref. [22J. 
Such a discrepancy could be of physical origin, due to the actual inequivalence of the 
0(3) and CP 1 models, contrary to the standard wisdom; or it could be of technical 
origin, due to shortcomings of the employed strategy in dealing with these models. 
The plan of the paper is the following. In Section [2] we briefly review the method of 
Refs. [201 EU E2J E3]. In Section [3] we describe the model of interest, discussing in 
particular the theoretical prediction for the critical behaviour of the model at 9 = n 
in the continuum, and working out the consequences for the observables relevant 
to our method. In Section @] we describe the 0(3)eNLaM on the lattice, and we 
discuss the results of our numerical simulations. Finally, Section [5] is devoted to our 
conclusions and to an outlook on open problems. Details of the numerical analysis 
are reported in the Appendix. 
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2 Theories with topological term in the action 
and the method of scaling transformations 

In this Section we briefly describe the relevant formalism and notation that will be 
used in the rest of the paper. The partition function of a theory with a topological 
term in the action is of the general form 



Z{9) = V(pe 



,-S[<p]+i6Q[<p] _ P -VF{8) 



where <fi denotes the degrees of freedom of the model, T><j) is the appropriate func- 
tional measure, S is the non-topological part of the action, and Q is the quantised 
topological charge, taking only integer values; moreover, F(9) is the free energy 
density and V the volume of the system. Clearly, Z{9) is a periodic function of 9, 
Z(9+2tt) = Z{9). In the interesting cases, the integration measure is invariant under 
parity (V), and S is P-even, while Q is V-odd. As a consequence, Z(—9) = Z(9); 
combining this with periodicity, we have that Z(ir + 6) = Z(n — 9). 

While at 9 ^ 0, 7r parity is explicitly broken, at 9 = n any "P-odd observable has 
vanishing expectation value in a finite volume; nevertheless, a phase transition may 
take place at this point. A convenient order parameter is given by the topological 
charge density, 

.(Q)„ 1 dlogZ(9) 8F{9) 
° {9) = ~ l — = ~V 89 = ~dT' (2) 
where we have introduced the notation 

{0[(j)]) ie = Z(9)- 1 [ D0e- 5M+i9QM 0[0] , (3) 



for the expectation value of the observable [</>]. In the limit of infinite volume, a 
nonzero value of 0(9 = it) indicates a first-order phase transition, while a divergent 
susceptibility O'{0 = n) indicates a second-order phase transition, and so on. 

In order to reconstruct the behaviour of the order parameter near 9 = ir using 
numerical simulations, one has to start from imaginary values of the vacuum angle 
9 = —ih, with /iGl. It has been suggested in Ref. [20] that a convenient observable 
is the quantity 

y(z) = -J^, z = cosh|, z>l. (4) 
V tanh | I 

It is immediate to see that under analytic continuation h — )■ i9 one has 

. (Q) ie 0(9) 9 

viz) = —i s = 5", z = cos — , Z<1. (5) 

yV ; V^tanf tanf ' 2 ' K J 

i.e., in terms of z the analytic continuation is simply an extrapolation from z > 1 to 
z < 1. Notice that y(l) = 2 ^^ fc , and y(0) = 0, with z = corresponding to 9 = 7r|]] 



1 One can have y(0) ^ only if the topological charge density diverges at 8 = ir, which seems 
unlikely. 
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The use of this observable is suggested by the antiferromagnetic one- dimensional 
Ising model, where the role of the (9-term is played by the coupling with an external 
imaginary magnetic field (for an even number of sites). This model is exactly solv- 
able, and one finds that y actually depends only on a specific combination of z and 
of the antiferromagnetic coupling F, namely y(z, F) = Y((e~ 4F — l)~^z). Although 
this property is exclusive of the one-dimensional Ising model, nevertheless one can 
expect that a similar smooth relation exists between y(z) and y\(z) = y{e^z) also 
in other models with #-term. The assumption usually made is that y(z) is a mono- 
tonically increasing function of z, vanishing only for z = (i.e., the order parameter 
does not vanish for 9 G (0, 7r)); this is indeed the case for the models where the exact 
solution is known. The quantity y x is then a monotonic function y\{y) of y, with 
the property that y\ = at y = 0, so that starting from the smallest values of y 
that can be obtained by numerical simulations at real h, one can therefore reliably 
extrapolate towards y — y\ — 0, i.e., in the region corresponding to real 9 = —ih. 
This is the advantage of this method, based on scaling transformation, over other 
approaches that involve an uncontrolled analytic continuation from imaginary val- 
ues of 9. Having reconstructed y\(y) in this region, one can then easily reconstruct 
the order parameter at real 9. Clearly, the closer one gets to y — 0, the better 
the extrapolation is expected to be: this method is then expected to work well in 
situations where the density of topological objects is small, such as asymptotically 
free models at weak coupling. 

If one is interested only in the critical behaviour at 9 = it, it is possible to de- 
termine the critical exponent without explicitly reconstructing the order parameter. 
Consider the effective exponent 7a (y), 

7AU/) = rlog • (6) 

A y 

Assuming a critical behaviour O oc z e near z = 0, i.e., O oc (n — 9) e near 9 = n, one 
immediately sees that y oc z e+1 near z = 0, and so 

7 = lim7 A (y) = -limlog — = l + e. (7) 

3/-V0 A z l+e 

Analogously, assuming that y\{y) is analytic at y — 0, one can obtain 7 from the 
relation 7 = f l O g(^| y=0 ). 

The method outlined above has been checked against explicitly solvable mod- 
els, and successfully applied to models where the exact solution is not known (see 
Refs. [20j [21j [22l [23]). One implicit assumption of this method is that the function 
y\(y) has a "reasonable" behaviour near y = 0, i.e., it can be well approximated by 
polynomials, or ratios of polynomials, or other "simple" functions. If this is the case, 
the critical exponent can then be obtained with fair accuracy. What has not been 
done yet is the evaluation of the impact of logarithmic corrections on the reliability 
of the extrapolation. The result Eq. ([7]) holds independently of logarithmic correc- 
tions to the critical behaviour, i.e., it holds even if O oc z e \og(l/ z)~P; nevertheless, 
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the way in which the limit is approached in this case can make the extrapolation 
more difficult. This issue will be discussed further on in the next Section. 



3 The 0(3) nonlinear sigma model with a topo- 
logical term 

In this Section we briefly recall the main properties of the 0(3) nonlinear sigma 
model with a topological term (0(3)gNLaM) in two dimensions, and we work out 
the consequences of the expected critical behaviour at 9 = tt for the method of 
scaling transformations described in the previous Section. 

3.1 Critical behaviour at 9 = tt 

The degrees of freedom of the 0(3)6iNLcrM in two dimensions are real three-compo- 
nent spin variables s(x) of modulus one, s(x) 2 = 1, "living" at the point x € M 2 . 
Expectation values are defined in terms of functional integrals as follows, 

(O[s\) i0 = Z(6)- 1 J Vse- s ^ +l6Q ^ 0[s\ , Z{9) = J Vse- s ® +i0Q ® , (8) 

where the measure is given by T>s = Y[ x d 3 s(x)5(l — s(x) 2 ). In the continuum, 

S[£\ = ^J d 2 x d^s(x) ■ d,s(x) , (9) 

and the topological charge Q[s\ is given by 

Q[^\ = ^ [d 2 x s(x) ■ e^d^x) A d u s(x) . (10) 

87T J 

Here [i, v — 1,2, and sum over repeated indices is understood; the antisymmetric 
symbol e^ v is defined as e 12 = — e 21 = 1, e 11 = e 22 = 0. 

While at 9 = the theory possesses a mass gap, it has been argued that as 
6 — > 7i the mass gap m(8) vanishes with the following behaviour 



2 

m(9) oc [tt — 6*| 3 



log 



7T 



9 



- — i 

8<iT,8~ir \ 7T — 9 I 



where we have neglected subleading termsH This prediction follows from the fol- 
lowing considerations for the continuum theory (see Refs. [121 122 I2S] ) . Near 9 = tt, 
the effective action for the 0(3) sigma model is given by the 377(2) i Wess-Zumino- 
Novikov-Witten (WZNW) model p31 [HJ [15] , with a marginally irrelevant, parity- 
preserving perturbation, and a relevant, parity-breaking perturbation, whose cou- 
pling g is a function of (tt — 9) that vanishes at 9 = tt . Renormalisation-group 



2 From now on we will work in the interval 9 £ [0, tt], so that we can discard the absolute values. 
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arguments relate as follows the coupling and the correlation length £ of the sys- 
tem El 



r ocdQoger 1 x [1 + ((log^r 1 )] ; 



9 



neglecting sub leading terms, as m oc £ one finds 



m oc 5 3 log - 
9 



(12) 



(13) 



It is usually assumed that g oc (tt — 9) + . . ., so that Eq. (TTIi) immediately follows. 

Following Kadanoff, one expects that near the critical point 9 = n the free energy 
density F{9) is proportional to the square of the inverse of the correlation length, 
that in turn is proportional to the inverse mass gap, so that 

HO) oc ^ oc m{9f . (14) 

The order parameter for parity breaking O(0), defined in Eq. (j2]), is therefore ex- 
pected to show the following behaviour near 9 = tt, 

-i 



0{9) oc —±L oc (tt 



log 



1 



TX-9 



(15) 



where we have neglected subleading terms. This behaviour is conveniently rewritten 
as follows in terms of the variable z = cos ■ 



0(9) oc z't log - 



1 



2 v*<l), 

-1 



Z 



z < 1 . 



(16) 



The critical behaviour is therefore a second-order phase transition, with recovery of 
parity, with a critical exponent e = |. 

For future utility, it is useful to work out the first correction to the leading 
behaviour Eq. ( !T3|) . This does not require the knowledge of the 0((log£) _1 ) terms 
in Eq. ( TT2|) ; we have 



m oc g$ 



log 4 

9 



1 



3 log log | 



+ 



8 log i logi 



f r (log | 



;i7) 



where the function r(x) is of the form r(x) = + (r± log(x) + r^) /x + . . ., and we 
have omitted subleading terms at large x. We will assume that subleading terms in 
the relation g oc (n — 9) + . . . are suppressed as powers of n — 9, so that they can be 
safely ignored in the analysis of the following subsections. Using Eq. ffl7|) . we find 
for the order parameter 

-i 



O oc z 3 I log - 

z 



3 lo g logl J_ 
4 logi log* [ gJ 



where again r(x) is of the form r(x) = r + (fi log(a;) + r 2 )/x + 
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3.2 Effect of logarithmic corrections on the effective expo- 
nent. Extension of the method 



We can now work out how the predicted behaviour of the order parameter near the 
critical point reflects on the effective exponent 7a defined in Eq. ([6]). Using Eq. ( TT8l) . 
one finds that near z = the quantity y(z) has the following behaviour: 



-l 

4 



y(z) = y zi ( log - ) c(z) , (19) 



z 



where yo is some constant and c(z) = 1 — | log lo f z + . . . (see Eq. (ITS]) ), where the 

° z 

dots stand for subleading terms. We have therefore for ^(2) 



y x {z) = e$*y z* f log^ - -J c(e^) 

-1 



and so 



e^(l + ^_+...)y(z) 



, , 4 2, / A 1 \ 41 
7A(y) = TT + T log 1 + -—T +... = - + — T + (21) 



3 A \ 2 log 1 7 3 log 



Taking the logarithm on both sides of Eq. (1191) . we find to leading order log - 
I log + . . ., and plugging this into Eq. ( 12~TT) we finally obtain 



^- 3 { 1 + ^) + °{^i)- < 22 » 

The derivative of this function is infinite at the origin: as a consequence, the effective 
exponent changes abruptly at very small y, going from 7a — | — 1.33 at y = to 
7a — 1.43 at y = 10~ 6 (see Fig. [T|). From a practical point of view, this behaviour 
makes it very hard to obtain the correct extrapolation from numerical data: one 
would need high precision dcitci cit very small values of y in order to figure out the 
logarithmic behaviour. 

It is possible to modify the method of scaling transformations discussed above, 
in order to reduce the effect of the logarithmic corrections. Indeed, it suffices to 
consider a new function y(z), obtained by multiplying y by an appropriate factor, 
designed to cancel the logarithmic corrections at 9 = n. A convenient choice is 



y(z) = viz) cosh — log ( 1 H r- j = — cosh — log f 1 H r- 

yy ! yy ! 2 & y coshfy Vtanhf 2 & y coshf 



(23) 
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1.52 

1.5 - 
1.48 - 
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1.34 




(a) 



(b) 



Figure 1: (a) Theoretical prediction for 7a (y) up to order O (\og ^yj- 
(b) Theoretical prediction for 7 a (y) up to order O ( ,^ \ v( fl 



It is immediate to show that the extra term behaves as 



cosi 4 iog 



—> cos - log 1 H 

2 \ COS 



->■ 2 log - ; 



(24) 



therefore, near z = we have that oc zi, without logarithmic correctionsjfl 
Compared to similar functions yielding the desired logarithmic term, this choice has 
the advantage that the behaviour Eq. (|24|) of the extra factor has only corrections 
of order 0(z 2 ) near z — 0, and that the large-h behaviour of y is the same as that 
of the topological charge density, so avoiding possible distortions in the numerical 
analysis. Moreover, the extra factor is a monotonically increasing function of z, so 
that it cannot modify the monotonicity properties of y(z) (which is assumed to be 
a monotonically increasing function for the whole method to work). 

It is straightforward now to work out the theoretical prediction for the behaviour 
of the new effective exponent 



Clearly, 



- f-\ - 2 i y x 

i\\y) = T-iog — ■ 

A y 



7 = lim7 A (?/) = 7 + l = e + 2 = - 

y^o 3 



(25) 



(26) 



More generally, if the order parameter behaves as O oc z e (log |) , one can define 



(Q)h 

Ftanh| 



in order to take care of logarithmic factors. One has that y(z,p) oc z e+1+p near z 
logarithmic corrections. 



0, without 
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Near z = 0, we have that (see Eq. (TI8]) ) 



3 log log - 
1 ^ 

4 log! 



(27) 



where y is some constant, and also 



V\{z) 



A 7 _ 7 



A 7 _ 7 

e 2 3 1/ 2 3 



i 3 log(log | - |j 
4 logi-f 



+ . . . 



1 



3 / log log | A log log | 
4\ l oe I 2 



log* 



+ 



(28) 



A 7 

e 2 s 



3 A log log - 

1 - TTTT + ■■■ 



4 2 (log 



z/0) 



where the neglected terms in the last passage are of order 0([log(l/,2)] -2 ). Taking 
logarithms on both sides of Eq. ( 1271) one immediately sees that to leading order 
log ~ = I log 4 + . . ., and so one finds that@ 



_ _ 7 _ 49 loglog | 

ix{y) -3-12 



+ 



(29) 



Although there still are logarithmic effects in the approach to the limit value, the 
"jump" of the function between y = and y = 10 -6 is half as much as that of 
7a predicted above (see Fig. [Q. Moreover, it is easy to see that corrections to 
the leading-order relation between log^ and log 4 are vanishing as y — > 0, i.e., 

log - z = flog 4 + o(l), while in the relation between log | and log^ there are also 
subleading but divergent terms as y — > 0. The bottom line is that the use of 7a and 
y instead of 7a and y is expected to improve the numerical analysis. 

In concluding this Section, we want to add a few remarks. First of all, we want 
to stress that the results of this Section are expected to hold in the continuum limit, 
and they are based on the fundamental assumption that the critical theory at 9 = tt 
is the 577(2) WZNW model at topological coupling k — 1. We notice also that the 
prediction for the quantity y(z) has been derived using the behaviour of the mass 
gap near the critical point 6 = tt, i.e., near z = 0, so that it is not expected to 
hold for z > 1, where numerical simulations are feasible. On the other hand, due 
to its expected smoothness, the prediction for y\(y) should hold more generally in 
the region of small y, that is accessible to numerical simulations at sufficiently small 
values of the coupling. 
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(a) (b) 

Figure 2: (a) A unit square of the direct lattice, i.e., a site of the dual lattice, 
(b) Spherical triangle corresponding to the spins s±, s*2 and S3. 



4 Numerical simulations on the lattice 

In this Section we describe the setup of our numerical simulations on the lattice, and 
we discuss our results on the critical behaviour of the two-dimensional 0(3)eNLaM 
at 9 = 7r at finite lattice spacing. 



4.1 The 0(3) e NLcrM on the lattice 

In order to compute numerically the functional integrals Eq. (jSJ), one replaces the 
continuum by a square lattice A of finite size V, properly discretising the action. 
The simplest choice for S is 

1 2 

S ^ ^-JLJLi 1 - s » ■ s > + A)] = W + PS latt [s\ , = l/g 2 . (30) 

" x£A 1-1=1 

Here jl is a unit lattice vector in direction \i. The lattice action <Si a tt[s] is identical to 
the energy of the Heisenberg statistical model, so that the resulting expression for 
Z{9 = 0) gives (up to an irrelevant constant) the partition function of this model at 
temperature 1//3 (in units of the Boltzmann constant). As regards the topological 
charge, we have used the geometrical definition of Ref. 



Q g eom[s] = q(x*) , q(x*) = — [(aA)(si,s 2 ,s 3 ) + (crA)(si,s 3 ,Si)), (31) 



47T 



where x* are sites of the dual lattice A* (i.e., squares of the direct lattice A), and 
Si = s(xi) are the spin variables living on the corners of the squares (ordered 



4 This result holds with a milder assumption on the relation between g and it — 9, namely that 
g oc 7r — 9 + . . . with subleading terms suppressed with respect to l °f^° S ■ 
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p 




statistics 


0.9 


100 2 


2 • 10 6 


1.2 


100 2 


2 • 10 6 


1.5 


100 2 


2 ■ 10 6 


1.6 


200 2 


4- 10 6 


1.7 


350 2 


2 ■ 10 6 



Table 1: Details of the simulations. 



counterclockwise starting from the bottom left corner, see Fig. [2] (a)). Here we have 
denoted by (aA)(si, S2, S3) the signed area of the spherical triangle having as vertices 
s*i, s 2 , and s 3 (see Fig. [2] (b)): the absolute value of the area A and its sign a, i.e., 
the orientation of the spherical triangle, are given respectively by 

A = «i + a 2 + a 3 - ir , o = sign [s\ • (s 2 A s 3 )] , (32) 

with <y,i the angles at the corners of the spherical triangle; the two terms q(x*) = 
qi(x*) + (^(x*) in Eq. (131 p correspond to the two triangles in which each square on 
the lattice is divided. In terms of the spin variables one has 

exp (|(crA)} = [1 + si ■ s 2 + s 2 ■ s 3 + s 3 • si + isi ■ (s 2 A s 3 )] , 

p 2 = 2(l + s 1 -s 2 )(l + s 2 -s 3 )(l + s 3 -s 1 ). ( j 

Except for the exceptional configurations 

s*i ■ (s* 2 A S3) = , 1 + si • s 2 + s 2 • s 3 + s 3 • si < , (34) 

for which the topological charge is not defined, one has a = ±1 and A < 2tt. One 
verifies directly that <3 ge om has the correct continuum limit; moreover, imposing 
periodic boundary condition, it necessarily takes only integer valuesjf] 

Regarding the numerical simulation of the system, the non-linear dependence of 
Qgeom on the spins makes it hard to envisage fast algorithms; we have therefore used a 
simple Metropolis algorithm, supplemented by a "partial overrelaxation" algorithm 
to accelerate the decorrelation between configuration. This "partial overrelaxation" 
algorithm simply consists in proposing the usual overrelaxation step used when 
simulating the model at 9 = to the Metropolis accept/reject step. This algorithm 
turns out to be rather efficient, especially when is large and the topological content 
of configurations changes rarely; notwithstanding its simplicity, it turns out also to 
be very effective in decorrelating configurations. 

5 It is worth mentioning that the Mermin- Wagner- Hohenberg theorem [27l [28l [29] , that forbids 
the possibility of spontaneous magnetisation in the model at 9 = 0, can be easily extended to 9 / 
if the geometric definition Q g0 om of the charge is used. 
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y 



Figure 3: Plot of the effective exponent 7^ (y) for (3 = 1.5, together with the result 
of a Bayesian fit at fixed 7 = 7a (0) (solid line, Tab. [2] (left)) and with free 
7 (long-dashed line, Tab. [2] (right)). 

4.2 Numerical analysis 

We have performed numerical simulations of the 0(3)eNLcrM at various values of 
the coupling. For each value of /3 we have chosen 45 values of h, in such a way that 
the topological charge was measured for both z = cosh | and z\ = ez z = cosh 4p we 
used A = 0.5. The (real) values of h = i0 that we used lie below the line of (possible) 
phase transitions determined in Ref. |30j, so that the region where our simulations 
were performed and the real-# axis belong to the same analyticity domain in the 
complex-6 1 plane. The statistical error on the topological charge has been determined 
through binning. The lattice size was chosen in order for the finite size effects to be 
negligible (see Tab. [[]). 

We have then analysed the results for the effective exponent 7a by means of 
Bayesian fits [31] , based on the theoretical prediction described in Section 13 .21 In 
a nutshell, a Bayesian fit takes into account our knowledge (the so-called priors) 
about the parameters that we are fitting. A detailed account of the analysis can be 
found in Appendix [Aj here we will mainly discuss the results. 

The fits were based on the following general form of 7a, 

7a0/) = 7 + F{\og y j ,{af}) , (35) 

with F(x) as x — >• 00, that can be derived from the expected critical behaviour 
at 9 = 7r (neglecting terms that vanish as power laws). The values of the parameters 

(k) 

yo and a] are not determined by the theoretical analysis, and have been fitted to 
the lattice data. 

A first analysis has been carried out by fixing 7 to the theoretical value, 7 = |, 
and fitting the other parameters, starting with y only and progressively adding 



13 



0.001 0.002 0.003 0.004 0.005 

V 

Figure 4: Plot of the effective exponent j\(y) for (3 = 1.6, together with the result 
of a Bayesian fit at fixed 7 = 7a (0) (solid line, Tab. |3] (left)) and with free 
7 (long-dashed line, Tab. [3] (right), and short-dashed line, Tab. |4j. 

terms, in order of relevance. We have then used the information obtained on y to 
tune the priors for a second fit, letting all the parameters free to vary. The results 
are reported in Tabs. El [3] and [5J for (3 = 1.5, = 1.6 and (3 = 1.7, respectively. 
Finally, at (3 = 1.6 we have also tried a fit using information on ajj 1 , obtained from 
the fit at fixed in order to set the corresponding priors: the results are reported 
in Tab. HJ The results of the fit with the largest number of parameters are shown in 
Figs. [31 HI and El for (3 = 1.5, (3 = 1.6 and (3 = 1.7, respectively. 

From the results of the analysis described above, we conclude that the lattice 
data are compatible, within the errors, with the critical behaviour predicted from 
the WZNW model, at f) = 1.5, (3 = 1.6 and (3 = 1.7. On the other hand, data at 
(3 = 0.9 and f3 = 1.2 led to bad-quality fits when fixing 7 to the theoretical value, and 
to a value of 7 considerably smaller than the theoretical prediction when allowed 
to float (~ 1.9 for (3 = 0.9, ~ 2 for (3 = 1.2). This shows that the WZNW-like 
critical behaviour does not hold at small f3, breaking down at some critical value 
yet to be determined, but does not allow us to draw any conclusion on the details 
of what happens as one lowers (3. The problem is that our analysis assumes a given 
logarithmic factor in the critical behaviour of the topological charge at 9 = tc, rather 
than obtaining it from the numerical data. A few attempts have shown that if we 
vary the exponent of the logarithmic factor in y, as described in footnote [31 we still 
obtain a good fit but the value for the critical exponent 7 resulting from the fit 
changes, too. For this reason, we have not attempted a more quantitative analysis 
at p = 0.9 and (3 = 1.2. 

Summarising, our results are compatible with one of the two following scenarios. 

1. There is a critical value (3 = (3 C , above which the critical behaviour of the 
0(3)6»NLcxM at 9 = tt is exactly the one predicted by the WZNW model at 
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Figure 5: Plot of the effective exponent ■j\(y) for (3 = 1.7, together with the result 
of a Bayesian fit at fixed 7 = 7a (0) (solid line, Tab. [5] (left)) and with free 
7 (long-dashed line, Tab. [5] (right)). 

topological coupling k = 1. 

2. The critical behaviour becomes exactly the one predicted by the WZNW model 
only at infinite (3, but for (3 large enough, (3 > (3 C , the difference is not appre- 
ciable numerically. 

As for what happens at small (3, there are various possibilities. As one expects the 
system to undergo a first-order phase transition at 9 = n at strong coupling, in the 
case of the first scenario above there may be a sharp change in the nature of the 
transition from first order to WZNW-like second order, with the order parameter 
vanishing at 9 = it as O ex (it — 9) e with e = ~. It is however also possible that 
the nature of the transition changes continuously, i.e., the critical exponent e varies 
from to I as (3 increases, either reaching | at some finite value of (3 = (3 C , or only 
asymptotically, as in the second scenario above. 

5 Conclusions 

In this paper we have studied the critical behaviour of the two-dimensional 0(3) 
nonlinear sigma model with (9-term (0(3)6iNLerM) at 9 = ir, by means of numerical 
simulations at imaginary 9. Using the method of Refs. [20j [211 (22J [23], it is possible 
in principle to reconstruct the behaviour of the topological charge density for real 9, 
and so investigate the issue of parity symmetry breaking at 9 = 71. The theoretical 
expectation is that parity symmetry is recovered at 9 = 71 through a second-order 
phase transition, the behaviour at the critical point being determined by the SU(2) 
Wess-Zumino-Novikov-Witten (WZNW) model [T3], [HI [15] at topological coupling 
k — 1. 
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Assuming that this is the case, one can show that the method of Refs. [201 1211 
122} 123] is unlikely to yield the correct critical exponent, as the large logarithmic 
violations to scaling at the critical point make it difficult to reconstruct the critical 
behaviour from the numerical data. Assuming that the logarithmic violations are 
known, it is however easy to modify the method in order to overcome this problem. 
We have then been able to show that our numerical results for sufficiently large (3, 
i.e., for sufficiently weak coupling, are compatible with the expected WZNW-like 
behaviour at 9 = it, in agreement with previous numerical investigations [2| [16]. 

Several issues remain open. Although the modified method allows to take care 
of logarithmic violations, it is necessary to know them in advance in order for it 
to work properly. Infact, an incorrect assumption on these logarithmic violations 
could not be detected from the numerical analysis, and so would lead to an incorrect 
evaluation of the critical exponent. The bottom line is that our modified method 
can be used to test a theoretical expectation on the critical behaviour of a model 
with results from numerical investigation, but would not lead to conclusive results 
if one rather tried to determine the critical behaviour from the numerical data. 

For this reason, we have not been able to determine the critical behaviour of 
the 0(3)eNLcxM at smaller values of /3, although we have been able to exclude 
that it is the same as in the WZNW model. Also, due to the numerical errors, we 
are not able to tell if the critical behaviour is exactly WZNW-like, starting from 
some critical value of (3, or if it becomes WZNW-like only asymptotically. Further 
investigations are therefore required, in order to unveil completely the phase diagram 
of the 0(3) e NLaM. 

Being in agreement with Refs. [21 [16], the results of this paper are obviously in 
disagreement with those obtained in Ref. [22] for the CP 1 model. There are basically 
two possibilities: either the 0(3) and CP 1 model are not equivalent, contrary to the 
standard wisdom; or they are equivalent, and the results obtained in Ref. [22J are 
affected by the numerical problems related to the logarithmic violations, discussed 
in Section 13.21 In order to settle this issue, a new analysis of the numerical data of 
Ref. [22] is required, along the lines developed in this paper, which will be discussed 
in a forthcoming publication. 
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A Bayesian analysis 

The basic idea behind constrained fits, also called Bayesian fits, is to use the available 
information on a given physical problem in order to improve the fits to the numerical 
data. We will not go into the details, that can be found in Ref. [31] and references 
therein: in this Appendix we will briefly describe the method in order to define the 
relevant notation and terminology. Then, after slightly extending the theoretical 
analysis of Section 13. 2^ we describe its application to the problem at hand. 

A.l Constrained fits 

Constrained fits are performed in practice by minimising a modified chi-squared, 
the augmented chi-squared Xa Ug > defined as xL g = X 2 + Xprior> where x 2 is the usual 
chi-squared, and where Xprior contains extra information that is used to constrain 
the fit. If, thanks to prior theoretical knowledge, one expects the parameters to 
be fitted, call them ai, a 2 , . . . , a n , to be close to the values a 1; a 2 , • • • , a n within the 
ranges <Ti, cr 2 , . . . , a n , then one sets 



This approach allows to add as many terms to the fitting function as desired, con- 
trary to what happens with standard fits. The goodness of the procedure is judged 

by the convergence of the errors on the various parameters as the number of terms 

x 2 

is increased, and by the value of * U1 ? , which should be of the order of or smaller 



than 1. If this is the case, the resulting error at the end of the procedure is expected 
to give a reasonable estimate of both the statistical and the systematic errors on the 
parameters of the fit. 

A. 2 Sub leading terms in the effective exponent 

In order to perform a Bayesian analysis of our numerical data, we need to go a 
few steps further in the derivation of the theoretical prediction for the relevant 
quantities. Ignoring corrections that contain powers of z, one can show that the 
effective exponent *f\(y) is of the following form, 



In order to determine all the coefficients Y- (A), one should know the subleading 
terms in the relation between the mass gap and the coupling g, and moreover all 
the proportionality constants relating g with n — 9, the free energy with the squared 




(36) 




(37) 
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mass gap, and so on. However, the coefficients Yj. [ (A) are under control, as they 
are not affected by the subleading terms and by the unknown constants (assuming 
of course that g has power, or power/log corrections only, beside the term linear 
in it — 6). Explicitly, : (A) = Y^ = (— 7/4) k . One can therefore resum the 
corresponding terms in Eq. (1371) . obtaining 



7 7 log log 



7A(y) = 3<i-4 



Ml 

y 



(logf) 2 + |loglogf logf 



oo fc-1 flog log ^y] 

iiEE^)^ ■ (») 

10 § y k=l j=Q (logf > 1 



Notice that r/ fc) (A) = (l) k Y} k \\), with the coefficients r/ fc) (A) expected to be of 
order C(10 _1 ), as they contain factors of positive powers of |, and also powers of ~ 
(recall that we used A = 0.5 in our analysis). 

We have then tried a Bayesian fit retaining the first few terms of Eq. f|38|) . namely 
using the following function, 

log logf 

ix{y) = 7 < i - v 



(logf) +77 log logf logf 

+ 7/1(0) + fMy) + ffM + tfM } ■ (39) 



The powers of 7 have been chosen so that the coefficients in the functions 



a W ai 2) +4 2) loglog% 

/i(y) = , / 2 (y) = , _ C 3 y 



logfj (logf 

2 



4 3) + af } log log f + af (log log f Y 

/s(y) = — _ 4 — > (40) 



to y 



af + a[ 4) log log f + (log log f ) * + af (log log f ' ''' 

TT5" 
logM 



/ 4 (y) 



are at most of order 1, and actually expected to be of order C(10 _1 ), as explained 
above. Here rj = 7/(7 — 1): this relation is easily found by substituting 7 to the 
value |, obtained for the WZNW model, in the theoretical analysis. 
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A. 3 Details of the numerical analysis 



As already explained in Section 14.21 a first analysis has been carried out by fixing 
7 to the theoretical value, 7 = §, and fitting the other parameters, starting with y 
only and progressively adding terms, in order of relevance. As regards the priors, we 

(k) 

assumed a Gaussian distribution for the coefficients a)- , with central value and 
a m 0.1, while for yo we have chosen central value and a ~ 2. We report in Tabs. |2] 
(left), [3] (left) and [5] (left) the results for y and for the Xa U g/ n ° data, for (3 = 1.5, 
(3 = 1.6 and (3 = 1.7, respectively. 

We have then used the information obtained on y to tune the priors for a second 
fit, letting all the parameters free to vary. The central value for yo was chosen close 
to the result obtained with the first fit, with o equal to the corresponding error. 
The results for 7 and yo are reported in Tabs. [2] (right), [3] (right) and [5] (right), for 
(3 = 1.5, (3 = 1.6 and (3 = 1.7, respectively. At (3 = 1.6 we have also tried a fit in 
which we have used information on obtained from the fit at fixed 7, to similarly 
set the corresponding priors: the results are reported in Tab. HI 
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n° par. 


log 2/o 


Xaug 

n° data 


n° par. 


7 


log 2/o 


Xaug 

n° data 


1 


-2.3078(22) 


12 


2 


2.3032(18) 


-2.051(17) 


2.1 


o 
Z 


-2.699(13) 


.64 


3 


2.3360(46) 


-2.726(48) 


.75 


3 


-2.685(31) 


.63 


4 


2.3339(87) 


-2.69(14) 


.75 


4 


-2.65(20) 


.63 


5 


2.3462(91) 


-2.42(15) 


.62 


5 


-2.55(25) 


.62 


6 


2.343(11) 


-2.38(16) 


.61 


6 


-2.55(26) 


.62 


7 


2.348(12) 


-2.36(18) 


.60 


7 


-2.44(28) 


.61 


8 


2.346(11) 


-2.25(18) 


.57 


8 


-2.39(32) 


.60 


9 


2.344(13) 


-2.23(20) 


.56 


9 


-2.39(32) 


.60 


10 


2.346(15) 


-2.24(21) 


.56 


10 


-2.36(34) 


.59 


11 


2.347(14) 


-2.21(20) 


.55 


11 


-2.31(34) 


.59 


12 


2.345(12) 


-2.19(21) 


.55 



Table 2: (Left) Result of a Bayesian fit at j3 — 1.5 with 7 fixed. 

(Right) Result of a Bayesian fit at f3 — 1.5 with 7 free. 



n° par. 


log 2/0 


x 2 

Aaug 

n° data 


n° par. 


7 


log 2/0 


y 2 

Aaug 

n° data 


1 


-1.3427(47) 


2.9 


2 


2.3222(20) 


-1.136(40) 


2.6 


2 


0.08(18) 


2.2 


3 


2.3633(61) 


-2.935(87) 


.62 


3 


0.13(32) 


2.2 


4 


2.367(10) 


-2.99(14) 


.60 


4 


-2.30(34) 


1.2 


5 


2.3711(86) 


-2.77(22) 


.55 


5 


-2.59(26) 


1.1 


6 


2.373(10) 


-2.79(18) 


.55 


6 


-2.96(40) 


1.1 


7 


2.375(10) 


-2.76(21) 


.54 


7 


-2.83(44) 


1.1 


8 


2.374(14) 


-2.71(28) 


.53 


8 


-2.69(37) 


.99 


9 


2.375(14) 


-2.72(28) 


.53 


9 


-2.80(37) 


.95 


10 


2.377(15) 


-2.73(29) 


.53 


10 


-2.82(41) 


.95 


11 


2.377(15) 


-2.72(29) 


.53 


11 


-2.69(36) 


.93 


12 


2.378(17) 


-2.73(32) 


.53 



Table 3: (Left) Result of a Bayesian fit at /3 — 1.6 with 7 fixed. 

(Right) Result of a Bayesian fit at j3 — 1.6 with 7 free. 
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n° par. 


7 


log yo 


Xaug 

ti data. 


2 


2.3050(15) 


-2.230(47) 


11 


3 


2.3627(61) 


-2.927(88) 


.58 


4 


2.3597(67) 


-2.892(91) 


.56 


5 


2.3601(71) 


-2.87(19) 


.56 


6 


2.3581(76) 


-2.83(19) 


.55 


7 


2.3586(91) 


-2.82(21) 


.55 


8 


2.3580(90) 


-2.79(23) 


.55 


9 


2.357(10) 


-2.77(25) 


.55 


10 


2.357(12) 


-2.78(25) 


.55 


11 


2.357(12) 


-2.77(26) 


.55 


12 


2.357(13) 


-2.76(28) 


.55 



Table 4: Result of a Bayesian fit at j3 — 1.6 with 7 free (2). 



n° par. 


log £0 


Xaug 

n° data 


n° par. 


7 


log^o 


Xaug 

n° data 


1 


0.200(13) 


7.3 


2 


2.3052(28) 


1.10(19) 


4.8 


2 


-2.740(52) 


1.3 


3 


2.3654(98) 


-3.35(15) 


.87 


3 


-2.63(10) 


1.2 


4 


2.367(11) 


-3.35(16) 


.85 


4 


-2.41(21) 


1.2 


5 


2.373(11) 


-3.07(19) 


.72 


5 


-2.39(24) 


1.2 


6 


2.373(11) 


-3.08(19) 


.72 


6 


-2.30(29) 


1.2 


7 


2.376(12) 


-2.99(20) 


.68 


7 


-2.21(29) 


1.2 


8 


2.375(11) 


-2.85(21) 


.63 


8 


-2.23(33) 


1.2 


9 


2.375(12) 


-2.86(21) 


.63 


9 


-2.18(36) 


1.2 


10 


2.377(12) 


-2.83(22) 


.62 


10 


-2.14(36) 


1.2 


11 


2.377(12) 


-2.77(22) 


.60 


11 


-2.12(38) 


1.1 


12 


2.375(12) 


-2.71(22) 


.58 



Table 5: (Left) Result of a Bayesian fit at /3 — 1.7 with 7 fixed. 

(Right) Result of a Bayesian fit at j3 — 1.7 with 7 free. 
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